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Abstract 


A  model  for  ionization  waves  in  a  D.C.  gas  discharge  is  developed  in 

■  i 

a  straightforward  manner  based  on  the  linearized  first  three  moment 
equations  for  positive  ions  and  electrons  and’  Poisson’s  equation.  Slab 

i 

symmetry  is  imposed.  The  predictions  obtained,  by  applying  this  model  to 
discharge  conditions  for  which  ionization  waves  have  beeh  observed  are 
in  good  qualitative  agreement  with  both  the  results  of  experiment  and 
the  predictions  of  other  theories.  The  effects  of  including  small  pertur-, 
bations  in  ion  temperature  and  electron  neutral  momentum  transfer  collision 
frequency  and  energy  transfer  collision  frequency  are  also  discussed. 
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LINEARIZED  THEORIES  OF  IONIZATION  WAVES 

i  ‘i 

i  >  : 

i 

:  I.  ;•  Introduction 

i 

l  J  .  1 

;  , 

The  positive  column  of  most  D.C.  gas  discharges  is  a  large  luminous 
>.  • 

Region  lying  between  the  anode  glow  at  the  anode  and  the  Faraday  dark 
space  toward  the  cathode,  (Fig.’  1).-  Often  for  discharge  conditions  of 
particular  interest,  because  of  their  applications  in  laser  technology, 
the  column  can  be  seen i at  any  particular  instant  to  have  a  regularly 

i 

striated  appearance ' (Fig .  1) .  When  firs.t  observed  these  waves  were  some- 
what  of  a  mystery;-, however,  in  the  past  fifteen  years  several  models  have 

been  proposed  to  describe  such  ionization  waves.  Among  these  are  the 

;  ! 

.models  of  Pekarek  (Ref  3) ,  Weissglas  and  Andersson  (Ref  14) ,  and  Swain 
and  Brown  (Ref  12) .  The  primary  basis  for  each  of  these  models  is  the 

1  •'  i 

use  of  a  too  species  three  moment  treatment  to  describe  small  variations 

i 

in  local  electric  fields,  number  densities,  drift  velocities,  and  temper¬ 
atures.  These  models  differ  in  the  exact  form  of  the  moment  equations 

I  j  > 

used  and  in  simplifying  assumptions  made.  This  report  explores  the 

i 

prediptions  of  the  straightforward  one  dimensional  too  species  three 
moment  model  of  discharge  disturbances  derived  in  Section  II  and,  in 

i 

Section  III,  compares  these  predictions  with  the  results  of  experiment 
and  with  the  predictions  of  the  theories  of  Pekarek  and  Swain  and  Brown. 
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II.  Theory 


The  Equilibrium  Positive  Column  (Ref  2:  238-251) 

This  discussion  of  the  equilibrium  positive  column  of  a  gas  discharge 
is  intended  to  establish  relationships  which  will  prove  useful  in  later 
sections.  The  axial  electric  field  strength  E  in  the  column  is  approxi¬ 
mately  constant.  A  direct  consequence  of  this,  seen  from  the  one  dimen¬ 
sional  Poisson  equation 


3E 

3x 


N  )  :  0 
e 


CD 


is  that  electron  and  ion  number  densities,  N-  and  N  respectively,  are 
approximately  equal  everywhere  along  the  longitudinal  coordinate  of  the 
column  x.  The  primary  charge  carrier  production  mechanism  in  the  column 
is  ionizing  collisions  between  fast  random  electrons  and  neutrals.  The 
primary  loss  mechanism  of  charge  carriers  is  ambipolar  flow  to  the  walls 
of  the  tube.  This  ambipolar  flow  is  set  up  when,  prior  to  equilibrium, 
electrons,  due  to  their  greater  mobility,  diffuse  to  die  tube  wall  much 
faster  than  positive  ions.  A  radial  field  then  exists  such  that  the 
wall  is  at  a  negative  potential  with  respect  to  the  rest  of  the  column. 
This  field  helps  further  induce  ambipolar  flow  by  tending  to  annul  orig¬ 
inal  differences  in  number  density  by  pulling  ions  toward  the  wall  and 
.repelling  electrons  (Ref  2:  143-145).  Current  in  the  column  is  also 

approximately  constant  with  most  of  the  current  being  carried  by  the  more 
mobile  electrons.  The  slower  positive  ions,  meanwhile,  balance  the 
electron  space  charge. 
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Since  current  is  constant  and  there  is  no  build  up  of  charge  in 
equilibrium,  ionization  rate  must  therefore  balance  ambipolar  loss  rate. 
Making  use  of  this  fact  Von  Engle  arrives  at  the  differential  equation 


d2N  A  1  dN  .  aN  n 

TJ  +  Far  r  0 


where  N  represents  either  electron  or  ion  number  density,  r  represents 
the  radial  coordinate  measured  from  the  longitudinal  axis  of  the  tube,  a 
represents  ionization  rate  per  electron,  and  DQ  represents  the  coefficient 
of  ambipolar  diffusion  (Ref  2:  144).  The  solution  to  Eq  (2)  is 


j^='jo(rv/S7Dp  (3) 
o 

where  N}  represents  the  equilibrium  number  density  along  the  longitudinal 
axis  and  Jq  represents  the  zero  order  Bessel  function.  Assuming  that  at 
the  tube  wall  fr  —  R'l .  Null  —  0*  En  f3l  is  satisfied  when 


g-  =  (c  /R)2 

a  1 


where  represents  the  first  zero  of  Jq  , 

Electron  kinetic  temperature  T  can  be  obtained  directly  from  Eq  (4) 
by  solving 


i7=-  v. : 2x"Vx 


where 
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3  1. 


a  =  Ap  /8q7iT  V. 

o  o  s  i 


-X 


(6) 


n  ~R>Te  ^Vi 

Da  ~  “q - T" 

o 


(7) 


and 


X  =  -A 


Vi 


(8) 


(see  Appendix  A,  Eq  (103)).  Here  A  represents  the  slope  of  the  ionization 

efficiency  curve  in  ion  pairs /m/Torr/volt,  p  represents  the  pressure  of  ' 

o 

the  neutral  gas  in  Torr,  u+  represents  the  ion  mobility  in  m2/volt/sec, 
q  represents  the  charge  of  an  electron  in  coulombs,  rn  represents  electron 
mass  in  kg,  and  represents  ionization  potential  of  a  neutral  in  volts. 
Once  T  is  found  from  Eqs  (5)  and  (8)  then  a  can  be  found  from  Eqs  (4) 

V 

and  (7). 

Equilibrium  charge  conservation  is  represented  for  the  column  by  Eq 
(2) .  Hie  equations 


(9) 


■  Ke(i)  -  V 


(10) 


express  equilibrium  conservation  of  axial  momentum  and  energy  respectively, 
for  electrons  (ions)  (Ref  4:  50,  61)  (Ref  2:  123).  In  Eq  (9) 

represents  electron  (ion)  drift  velocity,  y*^  represents  electron  (ion) 
mobility,  and  E  represents  axial  electric  field  strength.  I.n  Eq  (10) 

Ke(i)  rePresents  the  average  electron  (ion)  energy  lost  to  a  neutral  per 
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4  *  ' 

collision,  <V  >  r.-,  represents  the  average  electron  ’ (ion),  random  velocity, 

*  i .  ,  l  ' 

^e(i)  rePresents  electron  (ion)  mean  free  path  between  collisions  i 

i  • 

with  neutrals,  and  T  and  T\  are  neutral  ‘and  ion  kinetic  temperatures 

1  •  ,  !  \ 
respectively.  For  reference  in  a  latfer  section.  Eqs  (9)  and  (10)  will 

appear  as  ,  .  (s  1 


iV-  ■  > 

-  U  =  E 

’  (ID 

.  V 

q  i  e  ■ 

0 

i  i 

1  1  . 

,  .  m.v. 

\ 

1  +,U.  =  E 

1  , 

1  ; 

(12) 

1 

~q  U  E  =  v  T  ! 

no  e,  eo  e  t 

i  1 

(13) 

>  ,  i 

) 

'  q  U.E  =  v.  (T.«  -  T  )  . , 

,  0  t  lo  1  0 

.  1 

i 

(14) 

\ 

!  1 

*  1 

is  the  electron  (ion)  momentum  transfer  collision  frequency 

l 

is  the  electron  (ion)  energy  transfer  collision  frequency. 

<#  v 

H  * 


T  is  typically  small  with  respect  ,to  Tg  in  Eq  (13)  and  has  j  therefore 
been  omitted.  . 

1  i  s 

i  • 

The  Basic  Physical  Mechanism  of  Striations  i 

» 

According  to  Pekarek,  the  basic  factors  in  the  production  of  moving 

i  1  j 

striations,  the  type  of  positive  column  ionization  waves  of  interest  in 

i  i  t  : 

this  report,  are  the  dependence  of  ionization  rate  a  on  local  electron 

1  .  !  1 
temperature ,  the  production  of  space  charges  and  hence  electric  fields 

i  j 

due  to  the  different  diffusion  rates  of  electrons  and  ions,  and  the 

changes,  in  local  electron  temperature  caused  by  space  charge  fields 

’  *  ■  !  , 
(Ref  9:  741)  (Ref  10:  893).  Assuming  the  above  mentioned  dependence 


i  i 
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of  a  on  Electron  temperature ,  the  sequential  interaction  of  the  above  to 
,  produce  moving  striatio'ns  is  summarized  by  Lee,'  Bletzinger,  and  Garscadden 

as  follows:  ■  :  1  .  »  '  1  1 

| 

*  *  *  • 

1  !  :  !  ' 

'  The  striations  are  cpnsidered  to  occur  at  a  sufficiently 
high  gas  pressure  so  that  ambipolar  diffusion  conditions  are  1 
operative.  The  mobility  of  the  electrons  is  much  larger  than 
that  of  the  ions  and, a  disturbance  in!  the  concentration  will, 
after  a  ihort  time,  produce  a  positive  space  charge  at  a  region 
where  the  concentration  is  increased.  This  space-charge  elec- 

•'  trie  field  will  cause  a  decrease  in  the  electron  temperature 
(giving  a  dark  region)  in  the  region  to  the  anode  side  of  the 
original  disturbance,  consequently 'reducing  the  ion  density.  .  , 

In  turn;  this, dark  region  produces  a  negative  space  charge 
and  an,  increased  Te  (forming  a,  bright  region)  and  thus  an 
increased  ne  td  its  anode  side.  ,This  assumption  is  in  accord 
with  the  experimental  profile  of  the  moving  striations 

1  (Ref  11:  381).  ! 

v 

|  .  !  ■ 

I  .  ■  , 

From  the  above  explanation  it  is  apparent  that  striation- phenomena 

I 

can  be  modeled  as  small  perturbations  of  tire  ionization  rate:  a,  the 

1  .  ,  ' 

charge  carrier  number  densities  N.  and  N  ,  Electron  kinetic  temperature 

■  "  i  e  , 

T  ,  and  electric  field  E.  Axial  variations  in  these  quantifies  were  the 

"  1  .  i  « 

basis  of  th'e  original  Pekarek  theory  (Ref  6:  452) .  In  a  very  recent. 

1  ' 

model  by  'Swain  and  Brown  based  on  the  first  two  ion  moment  equations  and 

1  ,  i 

i  the  first  three  electron 'moment .equations,  small i axial  perturbations  of  , 

t  •  i  «  ‘  1 

electron  drift  velocity  Ue,  ion  drift  velocity  lb,  and^  vQq  with  respect 

i  .  , 

to  T  were  also  considered.  Swain  and  Brora  further  consider  N.  always 
e  ■  ■  t  i 

'equal  to  Ng,  assume  1^  to  be  negligible,  and  allow  radial  variation  in 

!  I  :  i 

their  basic  and  equilibrium  equations.  , They,. however,  consider  no  vari-  * 

l  *  i 

iation  in  radial  velocities  and  only  wave  like  axial  variation  of  the 

i 

remaining  quantities  in  their  first  order  linearized  equations.  They 

'  I  <  \\  ! 

make  no  attempt  to  perturb  the  ambipolar  loss  term’/'  (Ref  12:  •'  1383-1386). 

r  ! 

The  model  discussed  in  the  following  sections  is  strictly  a  one  dimensional 

*  »  , 

1  •  ,  1  »  » 

I  * 

,  1 
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model.  In  addition,  N.'  is  not  considered  equal  to  N  ,  an  attempt  is  made 

to  vary  the  ambipolar  loss  team,  and  perturbations  of  T^,  v_  with  respect 

to  E,  and  v  with  respect  to  E  are  considered. 

Go 

The  Three ‘Moment  Equations 1 

The  discharge  conditions  and  plasma  oscillations  discussed  in  this 
paper  are  limited  to  those  for  which  the  motion  of  positive  ions  and 
electrons  can  be  adequately  described  by  the  first  three  moments  of  the 
Boltzmann  equation.  Assuming  f1  the  distribution  function  for  the  ith 
species  to  be  isotropic  in  a  system  drifting  with  an  average  velocity  q, 
the  first  three  moment  equations  for  any  species  are 


-  *1,  j_  dv  VJt(H]coll.  C17) 


where  represents  the  position  space  coordinate  in  the  £th  direction, 
t  represents  time,  <F>^  represents  the  sum  of  external  forces  averaged 

^le  treatr..ont  discussed  in  the  remainder  of  this  section  closely 
parallels  that  of  D.  A.  Lee  (Ref  8:  1-19) . 
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over  velocity  space,  a0 .  represents  the  kinetic  stress  tensor,  m  repre- 
sents  the  particle  mass  of  the  species,  p  represents  the  pressure  measured 
by  an  observer  drifting  with  velocity  q^, 

=  E  j-  d7  v*f  (18> 

;  _ 

and  n  =  _  dv  f  (19) 

•'v 

Here  tlie  usual  repeated  subscript  notation  is  used  to  indicate  simulation. 
Eq  (16)  describes  momentum  transfer  in  the  £th  direction.  Assuming  f  to 

be  Maxwellian,  the  kinetic  stress  tensor'  a0 .  becomes 

> 

aZj  =  p6S,j  =  nf<S)lj  t20-5 

where  T  is  the  kinetic  temperature  of  the  species  (Ref  1:  114-121).  If 

there  is  no  applied  magnetic  field,  the  self  field  is  neglected  with 

respect  to  the  applied  electric  field  in  the  Lorentz  force  term,  and  all 

other  forces  are  neglected,  then  <F>^  reduces  to  zq^E^  where  z  represents 

charge  number  and  sign,  q  represents  basic  electronic  charge  and  E0 

o  « 

represents  the  &th  component  of  electric  field  strength.  With  these 


assumptions  and  letting 


q  =  L  <^(H)coii. 

4  V  / 


=  t  *  vt(l] 


coll. 


^  =  2  L  dv  Y74at.  coil. 
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o 


*  l 


*  v 


Eqs  (15),  (16),  and  (17)  can  be  written  with  slight  simplification  as: 


3n 

9t 


+ 


(24) 


3q*  +  3q* 
at  qj  ax. 
J  3 


3T  3T  2 

at  ax^  "  "5 


1  3 

nT 

"n  3x!. 

Lm 

“Ma 


Zq0E& 

m 


‘  !  P  +  n(qA 


(25) 


(26) 


Moving  striations  propagate  longitudinally  in  the  positive  column. 

In  order  to  simplify  Eqs  (24),  (25),  and  (26)  variations  in  q,  E  and  x 
will  be  considered  to  talce  place  only  in  the  longitudinal  or  x  direction 
where  x  increases  from  anode  to  cathode.  This  is  not  of  course  a  good 
assumption  in  terms  of  what  is  actually  happening  in  the  column  since 
.radial  variations  in  charge  carrier  velocity  and  electric  field  strength 
exist  and  undoubtedly  have  their  effect  on  local  electron  temperature 
and  particle  production  and  loss.  However,  the  solution  of  the  much 
simpler  one  dimensional  equations  does  provide  qualitative  results  and 
an  insight  into  the  more  complex  three  dimensional  problem.  Thus,  letting 
q  =  {u,0,0},  E  =  (E,0, 0} ,  and  x  =  {x,0,0}  Eqs  (24),  (25),  and  (26)  become 


3N  3U  3N 

6  +  N  ~  +  U  7; — — •  =  Q 
•  e  3x  «  y 


3t 


e  3x 


(27) 


3N.  3U.  3N. 

Fr+Ni3F+ui3F=Q 


(28)' 


3U 


3U_ 


„  —  T.  3N  .  3T  q  E  QU 

£  +  n  e  +  e  +  1 _ £  a  JlSL.  +  a  +  r~£ 

3t  e  3x  N  m  3x  m  3x  in  e  N 
e  e  e  e  e 


(29) 
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M  -,V 

X  * 


3U.  au.  '  T.  3N.  -j  9T.  Zq  E  QU. 

1  .  |  t  X  ,  X  1-1  1  _  Q  a.  A  4.  X 

at  ui  3x  N.m.  3x  m.  3x  ni.  Ai  N. 


8Te  .  2  .  sue  2 


sf  +  Ue  IT  +  t  To  fcf  -  -f  VcUe  *!Pe*  <3^  *  Te  ™ 


lU  2 
e  e 


.  3T.  7  3U.  ~  7  .... — 

1  +  +  4  T;  =  -4  +  4  P4  +  Q  -M-  -  T.  (32) 

1  3x  3l3x  3lll3l^3  lj 


(m.U. 2 


Eqs  (27)-  (32)  along  with  Poisson's  equation 


3E 

3x 


(33) 


are  a  system  of  first  order  non-linear  differential  equations  relating 
N  ,  N . ,  U  ,  U. ,  T ,  T. ,  and  E.  The  subscripts  i  and  e  denote  ion  and 

6  1  6  1  6  X 

electron  quantities  respectively. 

The  Collision  Terms 

Q  represents  particle  production  and  loss  due  to  collisions.  As 
previously  mentioned,  the  primary  charge  carrier  production  mechanism  is 
ionizing  collisions  between  fast  electrons  and  neutrals.  This  is  just 
simply  the  rate  of  ionization  per  electron  times  electron  number  density 
or  aN .  Since  the  primary  charge  carrier  loss  mechanism,  ambipolar  flow 
to  the  walls,  is  a  radial  phenomenon,  no  realistic  charge  loss  model  is 
available.  In  order  to  balance  charge  production  and  charge  loss  in 
equilibrium  a  contrived  loss  term  --(N.  +  N  ) ,  similar  to  that  used  by 
Pekarek,  is  included  in  this  analysis  (Ref  3:  857).  Here  x  represents 

the  mean  life  time  of  the  charge  carriers,  x  probably  depends  upon  both 
a  and  T  ,  however,  since  the  nature  of  this  dependence  is  unknown,  it 

v 

will  be  considered  constant  throughout  this  treatment.  Hence, 
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Q  -  +  V 


In  equilibrium  the  acceleration  of  the  electrons  (ions)  by  the  applied 
electric  field  must  be  equal  to  the  rate  of  their  momentum  transfer  to 
neutrals  during  collisions.  Hence  from  Eqs  (29)  and  (30)  • 


A  l-L 
e  m_ 


m  Zq  E 
A  1  -  0 

Ai  nu 


Therefore  from  Eqs  (11)  and  (12) 


A  =  -v  U 
e  -  e 


Ai  ■  -v+ui 


In  equilibrium  Eqs  (51)  and  (32)  combined  with  Eqs  (35)  and  (36)  become 


q  ®e  =  meAeUe  =  Pe 


-Zq  EU.  =  m.A.U.  =  P. 
1  111  1 


Therefore  from  Eqs  (13)  and  (14) 


P  =  -v  T 
e  eo  e 


P.  =  -v.  (T.  -  T  ) 
i  io  i  o 
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The  Characteristic  Form  of  the  Equations 

With  Die  collision  terms  modeled  as  in  the  previous  section,  Eqs 
(27)  and  (33)  are  a  system  of  quasi- linear  equations  having  the  form 


A.  •  +  B. .  =  F . 

13  dX  13  dt  1 


(43) 


where  A. .,  B. .  and  F.  are  functions  of  the  variables  y..  The  equations 
of  such  a  system  can  be  written  as  an  equivalent  system  of  equations 
each  of  which  involves  differentiation  in  only  one  direction  in  the  x-t 
plane  at  each  point  (x,t),  provided  there  exist  j  vectors  V  =  V-,  i  = 

1,  ... }  j,  such  that  j  linear  combinations  of  system  (42)  of  the  form 


V.A.  2. 
111  ax 


3yi 

+  «  V.F. 


(44) 


exist  where 


V,A. .  = 

1 


xviBii 


(45) 


for  real  X.  If  j  such  linear  combinations  do  exist,  then  these  equations 
can  be  used  to  form  an  equivalent  system  of  equations  each  having  the 
characteristic  form 


3yi  - 1 

X  ■  ■  J-  +  _ i 

ax  at 


Vij  ■  viFi 


(46) 


where  the  differentiation  is  in  one  direction  only,  along  the  curve 
(x  -  xq)  —  X (t  -  t  )  (Ref  5:  103-107)  (Ref  8:  10). 

The  XJs  associated  with  the  system  (27) -(32)  are  U  ,  U  ±  a  ,  U.  and 

_  6  6  X 

Uj  *  where  ag  =  /jF/iiT  and  a^  =  /-|  T^/nu .  Eqs  (27) -(33)  written  in 
characteristic  form  are:  • 
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2  3N  3T  f  *  m 

ITe?r+Ne3i:=-§TeS  +  !Pe-meAeUe^r§“VNe  '  ™ 

c  e  e  ■ 


2  9N*  u*.  r  o  Jll-  n 

5/ Ti  sr +  Ni  357  =  Ti<!  * 1  pi  ■-  Viui +  T  Sr  ■ V  Ni  C48) 

1  1  I  I  1  I 


T©  3Ne  9U  N  3T.  _  q  B  QU J  ,  ,  .  '  N  P 

me  *  Neae  aX±  +  HT  ax±  =  +Neae  in~  ‘  Ae  +  N~  "3  NeAeUe,  +  3.  “in” 


e  e 


e  e  e 


,  +  §  U  • 2  (49)  and  (50) 


Ti  9Ni  9Ui  Ni  8Ti  "  fa  E  QU.'I  .  o  N.P. 

5[S5f  4  Vi  m*  JT  3W  =  fVilf  +  *i  -  |  NiAi®i  + 

Aj-  111  II  1  *  1 


:  +  §  U.2  *  (S1)  (52) 

i 

W"ere  =  It  +  Uo(i)lx  md  S%7  "4  +  <».to  *.Wfc  .**“« 

Eq  (33)  is  in  characteristic  form.  1  1 


The  Simplified  Electron  Equations  s 

•  1  i 

In  the  electron  equations  the  time  derivatives'  may  be  neglected  for 
any  electron  variable  Y  since  for  phenomena  of  interest  ’ 


3Yq 

9Y ! 

6 

3t 

«•  U  ~  « 
e  3x 

o>.!±  ae>sr 

(Ref  6:  454).  This  is  equivalent  to  saying  that  both  , the  electron  drift 
velocity  Ue  which  is  usually  >  103m/sec  and  the  electron  acoustic  speed 

*  1  »  .i 

ae  usually  >  10sm/sec  are  much  greater  than  the  phase  velocity  of  the 
ioni Action  waves  which  is  usually  <  102m/sec.  Making  use  of  statement  , 
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(53)  and  using  the  'fact  that  U  2  «  a  2 ,  Eqs  (47) ,  (49) ,  and  (50)  can  be 

!  .  i  C  >  ,C 

reduced  to 


3N  '  v  N  v  UN  , 

ST  =  S  TT~  -me  -T—  '  T  \me 
.  e  e 


+  «Ne  -  IfNj  +  Ne)  i- 


8Ue ,  v-U/m»  '  2  3  q®.  ’4  1  U  > 

5x-  =  “ 17  “  '  I  V  t  |  ■ +-*♦  jt«Ne  -  i(N.  ♦  iy 


3T  1,  v  T  -  .  T 

«rv!-r^KE-  [“vltvyio- 

v  C  C 


1  > 


•  , 

where  the  relationships  expressed  by  Eqs  (54) ,  (37) ,  and  (41)  have  been 
'  substituted;  for  Q,  A  ,  and  P  respectively. 

I  V  C 


The  Nondimens iona'l  Linearized  Electron1  Equations 

In  A  normal  equilibrium  discharge  N  ,  N. ,  U  ,  U. ,  T  ,  T.  and  E  are 

'  C  1  6  1  6  1 

approximately  constant.  Thus  the  equilibrium  conditions  resulting  from 
. 1  Eqs  (33) ,  (54) ,  (55),,  and  (56)  are 


N.  =  N  =  N° 
i  e 


2. 

a  =  7 


TJ  • 

=  -q  E°  =  v  U  V 
U  o  o  •  -  e  e 

e 


i 


i 
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where  the  zero  superscript  denotes  equilibrium  value.  Since  no  radial 
variation  in  number  density  is  allowed,  N°,  the  average  nunber  density 
in  ary  cross  section  of  the  column,  is  taken  for  the  equilibrium  number 
density  for  both  electrons  and  ions  throughout  the  column.  The  variables 
N. ,  N  ,  U  ,  T  and  F,  are  assumed  to  have  the  following  form:  Y  = 

Y°  +  yex pfifex  +  6t)  where  Y  represents  ary  one  of  the  above  mentioned 
variables,  Yc  its  equilibrium  value,  and  t/exp(ikx  +  a  small  pertur¬ 
bation  of  Y  away  from  the  equilibrium  value  with  variation  in  space  and 
time  as  indicated. 

As  previously  mentioned  a  may  be  considered  to  vary  with  T  ,  v_  with 

v  ” 

E,  and  ve  with  both  Te  and  E;  therefore,  small  perturbations  in  these 
appear  here  as 

a  =  a  +  a*  £t  (60) 

v  =  v  +  v  o.  '  (61) 


v  =  v  +  v  <l  +  v'  tz 
©0  ©0  ©0  ©0 


8  9 

where  '  (prime)  denotes  T  _T  o  and  •  (dot)  denotes  E=Eo.  a'  can 

e  e~e 

be  confuted  directly  from  Eq  (6).  Swain  and  Brown  (Ref  12:  1383-1386) 

consider  v_  constant;  however,  according  to  Von  Engle 


y  «  E 


when  equilibrium  discharge  conditions  are  such  that  collisions  between 
electrons  and  neutrals  may  be  considered  elastic  (Ref  2:  124) .  This 

variation  can  be  used  provided  the  time  to  reach  microscopic  equilibrium 
(approximately  10" 6  sec)  is  short  compared  with  the  period  of  oscillation 


IJ 
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t  : 
\  / 


*  'V, 


of  striations  (approximately  10" 3  sec).  For  elastic  collisions  based  on 
statement  (63) 


(64) 


For  inelastic  collisions  v_  will  be  approximated  as  constant;  hence 
v_  =  0.  Plots  available  in  the  literature  similar  to  that  in  Fig.  2  are 
used  to  determine  which  type  of  variation  is  appropriate  for  a  particular 
gas  and  discharge  condition  (see  Table  I).  Based  on  Fqs  (10)  and  (13) 


From  Eq  (65) 


and 


(66) 

(67) 


Electron  mean  free  path  X  is  approximately  constant  for  discharge  con- 
ditions  of  interest  here  (Ref  2:  33,  34).  For  elastic  eler  ron-neutral 

collisions,  the  fraction  of  energy  transferred  k  is  constant  and  -  0. 
Plots  available  in  the  literature  similar  to  that  in  Fig.  3  are  used  to 
•determine  which  type  variation  is  appropriate  for  particular  gas  and 
discharge  conditions  (see  Table  I) .  ' 

With  the  above  assumptions  the  linearized  nondimens ional  electron 
equations  are 
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Figure  2.  An  Example’  of  Equilibrium  Electron  Drift 
Velocity  |Ue°|  as  a  Function  of  Reduced 
Field  E°/p°  for  an  Arbitrary  Gas. 


Figure  3.  An  Example  of  Equilibrium  Average 
Fraction  of  Energy  Transferred  per 
Electron-Neutral  Collision  k  as  a 
Function  of  Reduced  Field  for  the  Same 
Arbitrary  Gas. 
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ikne  =  [|  a  +  «;*]  te  -  l  au0  +  ^  ae  -  i*(ni  -  n  )  (68) 

^  e 

p  2 

ikue  =  -f  a  *  ?  “e*  ve)  ze  *  I  aue  '  TO  ae  '  T  ?~Cni  '  V  ^ 
v  e 

ikte  *  ('I  a  -  “e*]  *e  *  !  aue  +  I  ae  +  F&i  '  V  ™ 

v  e 

for  disdiarge  conditions  such  that  momentum  and  energy  are  transferred 

t 

through  elastic  collisions  and 

ikne  =  (I  a  +  “e*K  '  I  aue  +  (I  a  +  I  aHe  '  ?Kni  '  V  C7» 
v  ■  v.  e 


iku  = 
e 


4  a  +  4-  a'*  u  2)  t  +  1  aa 

5  5  e  e  e  5  e 


(3  A  2  •«) 
k-  a  +  ^  ate*  e 


-  y  F"Cni  '  V 

e 


ikte  -  ('I  a  -  V  !  aue  +  (f  a  ‘  !  “H0  +  TO<ni  "  V  (73) 

'  v.  e 


for  disdiarge  conditions  sudi  that  momentum  and  energy  are  transferred 
through  inelastic  collisions.  The  nondimens ional  quantities,  except  a,  in 

v  ne  «i 

the  above  are  the  variables  x,  k,  n  .  n. ,  and  t  denoting  feR,  — , 
at  e  i  e  R  N°  N° 

- ,  and - respectively  and  the  nondimens ional  constants  a'*,  t*,  k*. 

»,  o  t  o  e  e  . 

Ue  e  T  °R  U  0  Eo  .  / — i - 1 

•and  pe  denoting  a*  — — ,  re  k,  andyUg0  me/Te°  respectively. 

^e 

The  constant  a  is  defined  by  the  equations 
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v  R  q  E°R  m  U  °R 

a  =  -£<>-  =  -  2? _ _  — y 

ii  0  i  t  0  ,  t  fl 

ie  i  e  '  'e 


The  R  appearing  in  the  above  quantities  denotes  tube  radius.  , 
The  Nondimens ionai  Linearized  Ion  Equations 

1  1  r  i  , 

In  addition  to  Eqs  (57)  mid  (58)=  the  ion  equilibrium  equations 

1  ! 

resulting  from  (33),  ’(48),  (51),  and  (52)  ar<=:  i 


•  v.  (T.°  -  T  )  • 
q  E°  =  v . u.- °m.  ='-2jl~L: - <L_ 

0  ’ll  y  0  » 

ii 


(74)’ 


(75) 


i  ,  •. 

Assume  the  variables  NL ,  Ng)  U.,  and  T.  also  vary  as  Y  =  Y°  +  {/expCifex  +  <S t) 

1  . 

as  in  the  case  of  the  electron  variables,  v.  and  v.  do  not  vhry  with  E 

since  y  and  the  fraction  of  energy  transferred, through  collisions  are 
+  ,  •  1 

constant  for  ions  (Ref  2:  113 j,  114)  •.  Any  variation  .of  with  respect  : 

to  is  small  and  can  be  neglected  based  on  the  following:  ‘ 

3v. 


-SP,  =  6[v.  (T.  -  T  )]  =  v.  t.  -  aT 

1  1q  1  0  1  o  1  c)K 


T.=T.0^i0  ‘  T0^i  = 

ii 


V.  +  y  V- 
lo  2  lo 


T 

1  _JL. 

(Ji  0 

I  1 


£.  ~  V.  t. 
1  lo  i 


(76) 


With  the  above  assumptions  t^e  linearized  nondimens ionai, ion  equations 


are: 


*  x 

A  ' 


-§(ik  +  s)nt  +  (ik 


^i  "  "  I  fti  +  3  bui  "  (s  "  pi2 
\  »  > 


ai\ 


1 

(f  '  I  yi2)wCni  "  V  '  (77) 

*>  -’l 
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i  * 


s 


s  + 

1 

ik 

'  5 

ni  4  3  Vi 

s  + 

fltlj 

ik 

u-  + 

s  •+ 

fi  tM 

ik 

“ij 

viJ 

i 

yi 

- 

- 

- 

- 

- 

ti  = 


±  L.  e  +■  fi  +  iJ 

yi  (3  yi 


bu.  -  4  ft. 

i  b  i 


5  2  i  5 

5'i'i'i 


“i^e  +  W(ni  -  ne5 


(78);  and  (79) 


U.1 


Hie  nondimens ional  variables  -t,  s , 'u. ,  and  t.  are  respectively  £-4- 

<  i  .  it.,  ■  K 


R  ui 


t.  .  , 

4-n-,  -i— ,  and  -^-1  Except  for  b,  the  nondimens  ional  constants  f,  a!*, 
U.°  -U.°  T.°  •  1  ‘  '  .  1  , 

1  ,x  1  ,  v.  R  ctRT  0  tU.°  ,  U.°  ' 

t|,  and  y^  in  tlie  above  equations  are  - — ,  -j~— ,  and  1 


respectively. 


U.°  ‘  U.° 

l  l 


/§T  i>i 


t  . 


V.m.U.  °R  q  E°R 
-  +  i  i  _  o 


!  .  T^° 


T  3  i 
1 


(80) 


The  linearized,  nondimensional -‘Poisson's  equation  is 


ikee  ,=  n.  -  n 
i  e 


(81) 


1  ’  e  E° 

where  s  is  — . 

Rq  N° 


The  Dispersion  Relation  ,  '  . 

i  !  '  !  ' 

Tlie  electron  equations  either  Eqs  (68),  (69)-,  and  (70)  or  (71),  (72), 

and  ,(73)  and  Eq  (81)  are;  solved  simultaneously  to  find  t  ,  and  e  in  tenfis 
<  •  e 

of  n.  such  that  1  ■  ' 

1  i 


*  v 
*%  * 


te  =  c,Ck)ni 


(82) 
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e  =  C  (k)n. 
2  1 


where  C  and  C  are  complex.  Eqs  (77),  (78),  (79),  and  (81)  can  be 
1  2 

reduced  to 


ss  +  £ 

3  s  3 


ik-  (5  5  2|UC  Sc  ikc|l 

S3  9  1  i  i  3  2  T|  jj  l 


+  I  bU;5  +  [-  |  f  -  ik  -  s  t.  =  0  (84) 


r  p  lkeC  r  A  ' 

ik  ■  §  V®!*0! +  i pi2  IT ni +  I ik  -  lb  ui 

1  *• 


+  |f+ik  +  sti  =  0  (85) 


-  5T  C2  +  5  ”i  “i*C,  '  -if4  ni+  !  V  *  !  +  y  ui 


+  f  \  *  0  (86) 

H1 


This  system  has  a  solution  only  for  those  combinations  of  s  and  k  which 
make  the  determinant  of  the  matrix  of  coefficients  zero  (Ref  7:  157), 

By  setting 


*  v 
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As 


*  V 


1  ik  -  l  ct!*C  '  +  i  C  ^ 
3  ill  3  2  x? 

5  -r. 
I  lk 

1  1 

a  dispersion  relation  of  form 

D(s,k)  =  s3  +  B(k)s2  +  C(k)s  +  D(k)  =  0  (88) 

\ 

is  obtained  where  B,  C,  and  D  are  complex. 

Striation-like  behavior  can  be  inferred  where  real  s(k)  has  a  posi¬ 
tive  maximum  (the  s(k)  at  which  disturbances  are  propagated  exponentially 
at  their  maximum  rate)  or  where  real  s (k)  has  a  relative  negative  maximum, 
which  would  indicate  only  slight  damping  (the  s(k)  at  which  disturbances 
are  the  least  damped  exponentially  relative  to  those  around  them  but  not 
so  heavily  damped  as  to  be  undetectable) .  The  phase  velocity  -  - 

and  group  velocity  -  as  well  as  frequency  -  -|j~-  and  wavelength 

1~-  at  the  point  of  the  striation-like  behavior  can  then  be  computed  from 
the  dispersion  relation  and  compared  with  experimental  results. 

The  Equilibrium  Data 

To  solve  the  dispersion  relation  Eq  (87) ,  it  is  necessary  to  'compute 
the  nondimens ional  constants  k *.  a,  b»  f,  a'*,  a!*,  x  *.  x- *,  p.  and  p„ 
from  the  equilibrium  values  E°,  p°,  N°,  U.°,  U  °,  T°,  T  °,  T. °,  \>_,  v  , 
veo’  Vio’  a>  T>  a'’  ^  ^ese  values»  all  which  are  in  mks  units 
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4  •v. 
1  » 


unless  otherwise  specified,  are  obtained  for  this  treatment  from  a  combin 

ation  of  experimental  data,  basic  theory,  and  the  equilibrium  equations 

themselves.  The  gas,  N°,  T°,  R,  and  p  are  known  or  assumed  based  on 

0 

typical  experimental  values.  Listed  in  Table  I  are  quantities  useful  in 
computing  gas  discharge  equilibrium  conditions  which  can  be-  found  in 
Von  Engle  (Ref  2) . 

Table  I 

Quantities  Useful  in  Computing  Gas  Discharge  Equilibrium  Conditions 


Quantity 

Symbol 

Units 

Source  (Ref  2) 

Ion  Mobility  at  1 

Torr  °C 

+ 

■V 

10" *m2 /sec/volt 

Table  4.1,  p  114 

Axial  Electric  Field 
Strength 

E° 

volt/m 

Fig.  125,  p  245; 
Fig.  127,  p  247 

Slope  of  Ionization 
Efficiency 

A 

ion  pairs/m/Torr/ 
volt/electron 

Table  3.7,  p  63 

Electron  drift  velo¬ 
city 

Ue 

103m/sec 

Fig.  61,  p  124 

Ionization  Potential 

Vi 

volts 

Table  3.6,  p  59 

Fraction  of  energy 
transferred/electron 
■neutral  collision 

K 

Fig.  63,  p  126 

• 

3k 

3E 

m/volt 

Fig.  63,  p  126 

Another  excellent  source 

for  basic  equilibrium  data  is 

Brown  (Ref  4) . 

Ihe  ion  mobility  at  1  Torr  and  0°C  y^+  adjusted  for  temperature  and 
pressure  according  to  the  equation 


y+  =  y0+  m  jr  .  C89) 
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-.v  can  be  used  to  obtain  v+  and  U. 0  from  Eqs  (9) . and  (12)  when  E°  is  known 

experimentally.  The  empirical  formula 

v.  =  (7  x  106)p  (90) 

io  o 

is  used  to  obtain  v.  .  With  v.  ,  T°,  U.°,  and  E  known  T.°  is  available 

io  io  1  1 

from  Eq  (14) . 

From  Von  Engle  U  0  can  be  found  knowing  E°  and  p°  (see  Table  I) . 

w 

The  electron  mobility  and  hence  v_  may  then  be  found  from  Eqs  (9)  and 

*  (11).  T  °,  a,  and  a'  can  be  obtained  from  Eqs  (5),  (6),  and  (8)  by 
making  use  of  Von  Engle  (see  Table  I)  to  obtain  A  and  V^.  With  v_,  E°, 

■j  U  0 ,  and  T  0  now  known,  v  is  easily  found  from  Eq  (13) .  From  Eq  (58) , 

*  6  C  Go 

*  2 

’  t  =  — .  See  Appendix  B  for, sample  calculation  of  equilibrium  data. 

,  Hi®.  of  Solution  _ , 

To  solve  the  dispersion  relation  Eq  (87)  k  is  assumed  to  be  positive 

real.  Starting  with  some  initial  value  of  k,  numerical  values  for 

and  C  in  Eqs  (82)  and  (83)  are  determined  by  solving  simultaneously 
2 

Poisson's  Eq  (81)  and  the  electron  equations  either  (68),  (69),  (70)  or 

(71) ,  (72) ,  and  (73)  depending  upon  whether  the  value  of  E°/p°  was  such 

that  inelastic  collisions  are  important.  Knowing  the  values  of  k,  C  , 

and  C  ,  numerical  values  for  each  element  of  the  determinant  in  Eq  (87) 

2 

are  found.  Evaluating  this  determinant  at  this  point,  values  for  B,  C, 
and  D  in  Eq  (88)  are  obtained.  Knowing  its  coefficients,  the  cubic 
dispersion  relation  is  then  solved  by  Newton's  Method  to  find  the  first 
of  its  three  complex  roots.  The  quadratic  formula  is  then  used  to  find 
the  other  two  roots,  k  is  then  incremented  to  a  new  value  and  the  above 

*  V 

, .  process  is  repeated,  k  is  most  conveniently  varied  so  that  the  results 
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are  easily  plotted  logarithmically.  By  using  electron  Eqs  (71),  (72), 
and  (73)  and  setting  <*  =  0  the  effect  of  not  allowing  v_  and  to 
vary  with  E  is  easily  evaluated. 
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III.  Results 


The  Roots  of  the  Dispersion  Relation 

The  real  and  imaginary  parts  of  the  solutions  to  the  dispersion 
relation  Eq  (87)  for  a  mercury-argon  discharge  at  3  Torr  (corresponding 
roughly  to  the  experiment  described  in  Lee  oX  at,  (Ref  11:  378))  are 

shorn  in  Figs.  4  and  5  respectively.  The  equilibrium  data  used  to  model 
this  3  Ton  mercury- in- argon  case  is  included  as  a  sample  calculation  in 
Appendix  B.  Roots  a  and  c  are  heavily  damped  over  the  entire  range  of 

k.  Root  b,  however,  exhibits  a  negative  maximum  at  approximately  k  =  3. 

The  narrow  band  of  frequencies  around  this  point  are  those  for  which 
disturbances  are  the  least  damped  and  therefore  the  most  likely  to  propa¬ 
gate.  If  the  damping  or  amplification  per  cycle,  2ir  >  is  less  than 

one,  a  wave  will  be  considered  to  be  neutrally  damped  since  for  such 
waves  there  is  a  high  possibility  that  damping  or  amplification  may 

"v 

appear  as  a  result  of  the  rather  crude  approximations  used  to  obtain  the 
equilibrium  conditions  of  the  discharge  modeled  and  not  as  a  result  of 
the  model  itself.  In  fact,  highly  believable  adjustments  to  the  starting 
conditions  for  the  cases  examined  here  will  actually  cause  the  model  to 
predict  2ir  ~y|y  much  less  than  one  or  actual  neutral  damping.  Never¬ 
theless,  the  damping  per  cycle  predicted  by  the  model  in  this  case  is 

l. 34,  which  even  by  the  above  criteria,-  indicates  the  presence  of  sig¬ 
nificant  damping.  The  striation  frequency  and  wavelength  at  the  point 
of  least  damping  are  8.7  khz  and  1.8  cm  respectively.  Striation  phase 

and  group  velocities,  v  and  v  ,  are  162  m/sec  and  -253  m/sec  respectively. 

Jr  o 

Hie  wave  predicted  is  a  backward  wave  since  the  phase  and  group  velocities 
are  in  opposite  directions.  The  experimental  values  of  v^  and  v^  are 
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SO  m/sec  and  -50  m/sec  respectively.  The  striations  in  the  experiment 

j 

are  neutrally  damped  and  occur  over  a  narrow  range  of  frequencies  (Ref 
11:  383).  While  backward  waves  over  a  narrow  range  of  frequencies  are 

predicted  by  the  model,  the  group  and  phase  velocities  predicted  are  of 
the  same  order  of  magnitude  but  not  approximately  of  equal  magnitude  as 
indicated  in  the  experiment.  Additionally,  the  damping  predicted  by  the 
model  does  not  show  up  in  the  experiment. 

Hie  Results  of  Other  Theories 

■;  Figs.  6  and  7  compare  the  real  and  imaginary  parts  respectively  of  • 

the  striation  roots  of  three  different  models  for  the  same  mercury-in- 
1  argon  case  depicted  in  Figs.  4  and  5.  Root  d  is  computed  exactly  as 

!  '  root  b  except  that  v_  is  held  constant  with  respect  to  E  rather  than 

;  #  considered  proportional  to  as  in  root  b.  Root  d  is  included  for  two^ 

j  v  '  reasons.  It  illustrates  the  effect  of  the  additional  assumption  v_  «  E^ 

.  when  compared  with  root  b  and  compares  favorably  in  the  region  of  stria¬ 

tion-  like  behavior  wit! a  root  e  computed  by  Swain  and  Brown's  model  (Ref 
12:  1383-1386) .  Root  f  is  computed  by  the  Pekarek  theory  as  detailed 

,  .  in  Lee  at  aZ .  (Ref  11:  381,  382). 

Roots  d  and  e  predict  neutrally  damped  waves  with  amplification  per 
cycle  of  .457  and  .828  respectively.  Root  f  predicts  exponentially  grow¬ 
ing  waves  with  an  amplification  per  cycle  of  2.79.  Striation  frequency 
and  wavelength  at  the  points  of  maximum  amplification  are  3 . 5  khz  and 
2,1  an  for  root  d,  3.58  khz  and  2.1  cm  for  root  e,  and  6.3  khz  and  1.34 
cm  for  root  f .  Hie  corresponding  group  and  phase  velocities  for  roots 
d,  e,  and  f  are  respectively  80.6  m/sec  and  -81  m/sec,  79  m/sec  and  -89 
in/sec,  and  84.5  m/sec  and  -35  m/sec.  The  sharp  peak  in  each  root 
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Figure  7.  The  Wave- like  Root  of  the  Dispersion  Relations 
(Iro(s)  vs  k)  for  Different  Theories 
(Mercuiy- in-Argon  at  3  Torr) 
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indicates  a  rather  narrow  range  of  frequencies  for  which  wave  propagation 
is  allowed.  The  results  predicted  by  the  Section  II  model  with  v_  f  v_(E) 
and  the  Swain  and  Brown  model  are  very  similar.  Both  agree  well  with 
experiment  predicting  neutrally  damped  backward  waves  over  a  n.  rrow 
frequency  band  with  group  and  phase  velocities  of  approximately  equal 
magnitude.  Hie  Pekarck  theory  predicts  growing  backward  waves  over  a 
narrow  frequency  band  with  phase  and  group  velocities  of  unequal  but  of 
the  same  order  of  magnitude. 

The  Characteristic  Behavi or  for  Large  k 

The  characteristic  behavior  of  the  equations  can  be  seen  in  the 
large  fe  limit  of  the  imaginary  parts  of  roots  a,  b,  and  |c  in  Fig.  4. 

Such  behavior  is  expected  since  for  variation  of  the  fora  expfifex  +  it) 

Eq  (46)  becomes 

(Aife  +  ^VjB^y.  =  V.F.  ‘  (PI) 

As  fe  -*  00  the  right  hand  side  of  Eq  (91)  becomes  negligible  .  th  respect 
to  tiie  left  hand  side;  hence 

(Aife  +  ^V.B^y..  :  0  '  (92) 

m 

(Aife  +  l)  ~  0  (93) 


(94) 


As  fe  gets  large  the  phase  velocities  for  roots  a,  b,  and  c  respectively 
approach  +U.°,  +U.°  -  Jt  T.  °/nw  ,  and  +U.°  +/|  T.  °/m. .  In  contrast,  root 

X  1  »  5  1  1  1  W  1  1 

e  in  Fig.  7,  based  on  Swain  and  Brown’s  model  which  assumes  =  0  and 
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Ni  =  Ne  throughout,  lias  a,v  in  the  large  k  limit  o£  +Ui°  -  Tf/nn.' 

Swain  and  Brown's  other  root  ;(not  pictured)  has  a  phase  velocity  of  1 
+V  +  J'j  as  k  Both  roots  b  and  c  initially  parallel  their 

corresponding  Swain  and  Bram  roots;,  however,  at  k  ~  30  the  difference 
in  N.  and  Ng  as  expressed  by  Poisson's  equation  becomes  large  enough  to 
separate  the  roots  of  the  two'  .iffer.  t,  models  toward  their  separate 

i 

asymptotic  paths.  i  '  . 

,  !  ' 

Additional  Comments  '  *  > 

-  | 

I 

Results  obtained  by  applying  the, Section  II  model  to  numerous  other  1 

*  i  l 

discharge  conditions  were  generally  similar  to  those  discussed  for  the 

mercury- in- argon  case  above.  In  particular  good  qualitative  agreement 

was  obtained  with  both  the  Pekarek  theory  predictions  and  the  results  of 

'  :  *  1 

the  experiment  described  for  the  50  Torr  neon  diffuse  discharge  case 

l  J 

described  in  an  article  by  .Garscadclen  and  Lee  (Ref  13:  578).  The  Section 

«  1  !  : 

II  model  did  not  predict  striation-like  behavior  for  the  35  Torr  argon  1 

diffuse  discharge  case  discussed  in  the  same  article.  However,  when  the 

•  .  *  «  ‘ 

fraction  of  energy  transferred  'in  electron  neutral  collisions  k  was  . 

.  i  | 

considered  constant  with  respect  to  E  or  -v  i  v  (E) ,  the  predictions  ' 

J  Go  Go 

of  the  model  were  in  much  better  agreement  with  both  the  Pekarek  theory  , 

!  1  \ 

and  experiment.  The  overall  effect  of  letting  v  .and  v  vary  with  E  is 

«  '  |  *  t  j  , 

to  increase  the  damping  in  the  root  containing  the  striation-like  behavior. 

•  V  * 

The  striation-like  behavior  in  the  majority  of  cases  other  than  the  one 
for  mercury- in- argon  discussed  hero  appeared  in  the  root  whose  large  k 
phase  velocity  approaches  U.°  rather  than  lh°  -  Jt  TP/nw  as  in  root  b, 

Fig.  4.  The  cause  of  such  behavior  and  whether  there  is  any  corresponding  . 

:  i 

physical  significance  lias  yet  to  he  resolved,  , 
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■  IV.,  Conclusions 

■  ■,  ■  \ 

Because  of  the  oversimplification  of  the  physics  involved  and  the 

small  amount  of  data  available  from  the  actual  experiments  modeled,  any 

i  '  J 

comparison  made  between  the-  results  predicted  by  any  of  the  models  dis- 

1  *  .  i 

cussed  herein  and  those  of  experiment  must  be  highly  qualitative.  In 

’  i  ’  ,  (  ,  i  ‘  , 

this  highly  qualitative  sense,  'the  predictions  of  all. the  models  discussed 

!  1  .  .»  1  ,  -  1  ■' 

in  thej  results  agree  moderately  well  with  experiment.  The  better  agree- 

.  ,  ‘  i  •  '  !  ‘  .  ! 

ment  with  experiment  obtained  by  the  Section  II  model  with  v_  f  v_(E) 

i  .  -  ,  -  ,  i  . 

(root'd)  and. by  Swain  aid.  Brown's  model  (root  e)  than  by  the  Section  II  1 

5  ‘  I 

model  with  v_  =  (E)  (root  b)  indicates  that  for  the  3  Torr  mercury- 

argon  discharge  y_  is  perhaps  better  approximated  as  constant  rather  than 

-i  1  ■  ^  >  ,i  1 

proportional' to  E  *.  Tending  to  confirm  the  above^equi librium  data  in 

Brown  (similar  data  for  an  arbitrary  gas  is  depicted  in  Fig.  2) funlike 

■  ;  '  I  I 

that  in  Von  Engle  (see  Table  I)  .indicates  that  this  discharge  condition 

,  _  '  1  i  _1 

is  not  quite  such  that  y_  can  be  considered  proportional  to  E  ?  and  that 

1  i  1  ' 

y  is  in  fact  almost  constant  with  respect  to  E  for  this  operating  region 

•  (  *  I  * 

(Ref  4:  55).  t  , 

1  As  mentioried  in  Section  III,  the  Section  II  model  produced  a  similar 

.  1  1  ' 
but  more  damped  behavior  when  v_  and  were  allowed  to  vary  with  E  than 

*  *  i  1 

when  they  were  held  constant.  Tins  seems  physically 'consistent  in  both 

I  v  x 

cases.  In  the  v_  case  the  effects  of  a  small  change  in  local  E. would  be 

1  ’  i  •  'i 

partially  offset  by  the  effects  of  a  corresponding  change  in  local  mobil¬ 
ity.  In  the  vgQ  case  the  effects  of  a  small  change  ,in  the  local  fraction 

;  ! 

of  energy  transferred  due  to  a  small  change  in  E  would  hasten  an  opposing 

1  i 

change  in  electron  temperature  and  hence  electric  field.  Tlie  disagreement, 
of  experimental  results  with  severe  damping  predicted  by  the  Section  II 


i 
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model  for  the  briefly  discussed  35  Torr  argon  case  was  probably  caused 
^  by  the  use  of  too  large  a  value  for  k  when  computing  the  dispersion 

relation.  The  fact  that  better  results  were  predicted  when  v  was  held 
constant  with  respect  to  E  indicates  that  k  cannot  be  reasonably  inferred 
from  equilibrium  changes  in  k  with  respect  to  E  (see  Table  I).  It  is 
hoped  that  using  a  more  reasonable  value  for  k  will,  in  the  future,  pro¬ 
vide  better  results  for  discharge  conditions  where  inelastic  electron 
neutral  collisions  are  important;  however,  for  the  purpose  of  this  report 
the  Section  II  model  containing  the  assumption  that  =  v  (E)  based  on 
the  above  method  of  obtaining  k  is  useful  only  in  inferring  the  possible 
effects  of  such  variation, on  ionization  waves. 

The  fact  that  all  models  examined  predict  striation-like  behavior 
in  region  from  k  =  3  to  k  =  5  for  the  3  Torr  mercury-in- argon  case  indi¬ 
cates  that  they  are  generally  consistent  with  each  other.  This  is  not 

\  * 

surprising  since  they  are  all  based  on  some  linearized  form  of  the  moment 
equations.  For  the  cases  examined  the  Section  II  model  with  v_  f  v_(E) 
and  v  f  v  (E)  and  the  Swain  and  Brown  model  predict  results  which  are 
in  the  best  agreement  with  experiment.  Solving  the  moment  equations  by 
■  hot  assuming  N.  =  N  and  T.  to  be  negligible  is  of  some  interest  from  a 
mathematical  point  of  view.  However ,  the  fact  that  Swain  and  Brown's 
model  which  includes  these  assumptions  and  the  Section  II  model  with 
v  i  v  (E)  and  v  f  v  (E)  which  does  not  include  these  assumptions 
predict  almost  identical  striation  behavior  (roots  d  and  e) ,  indicates 
that  the  more  complicated  cubic  dispersion  relation  Eq  (87)  obtained  in 
the  case  of  Section  II  model  offers  little  advantage  over  the  simpler 
quadratic  dispersion  relation  obtained  by  Swain  and  Brown. 

4  v 
*  * 
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Appendix  A 


The  Ionization  Frequency  Equation 


Beginning  with  Von  Engle’s  equation 


“"irf  h1 

a  '  p. 


e 


£i)fcN(e)de 


where 


N(e)d(E)  = 


1  _£ 

2N  t  e  r  > 
— £_  e  md  — 
/fr  em  em 


f .  =  \~/K 

c  Jme  e 


ei  =  Vi 


em  =  Te 


Apn 

c  =  r-2-  Xa 
e 


C100) 


and  e  is  energy  (Ref  2:  293),  the  ionization  frequency  a  becomes 


a  =  Ap 


/8/m  ir  e  ^  ( 

»  e  m  J 


(e  -  e^)ee  de 


(101) 


Letting  e  =  +  X,  Eq  (100)  reduces  to 
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where  is  the  ionization  potential  in  volts  and  A  is  the  slope  of  the- 
ionization  efficiency  curve  in  units  of  ion  pairs/m/Torr/volt. 
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’  Appendix  B 

Sample  Calculation  of  the  Discharge  Equilibrium  Conditions 
for  the  3  Torr  Mercury- in-Argon  Case  Discussed  in  Section  III 

Values  for  N°  and  T°  are  assumed  based  on  typical  discharge  conditions 
to  be 

N°  =  1016  m"3 
T°  =  5.52  x  10’21  (400°K) 

R,  p  ,  E°  and  T  0  from  the  experiment  modeled  (Ref  11:  383)  are 

o  e 

\ 

R  =  .01  m 

p  =  3  Torr 
ro 

E°  ^  180  volt/m 
Te°  =  2.08  x  10"19  joules 

The  mercury- in- argon  mobility  y  +  at  300°K  and  760  Torr  taken  from  Brown 

0 

(Ref  4:  77)  is 

y  +  =  .00018  m2 /volt/sec 
0 

y  +  corrected  to  400°  and  3  Torr  is 
0 

y+  =  .0605  m2 /volt/sec 

From  Eqs  (9)  and  (11) 
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v,  = 


=  —o—  =  7.7  x  10 6  sec"1 


+  mHg’,+ 


From  Eq  (9) 


lb0  =  u+E°  =  10,8  m/sec 


Frcm  Eq  (90) 


v.  =  2,1  x  107  $ec-1 

io 


From  Eq  (75) 


'  q  E°U.° 

T.°  =  -2 - —  +  T 

Vio 


Ti°  =  5,53  x  10~21  joule;; 
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From  Von  Engle  (see  Table  I)  for  electrons  in  argon 


and  for  mercury 


U  0  =  -4  x  103m/sec 

v? 


V^  =  10.4  volts 


From  Eq  (11) 


EY 


v  = 


mil 
e  e 


=  7.9  x  io’9  sec"1 


From  Eq  (13) 
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From  Hqs  (4) ,  (5) , 


a 


Differentiating  Eq 


From  Eq  (58) 
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